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One goal in the quantum-walk research is the exploitation of the intrinsic quantum 
nature of multiple walkers, in order to achieve the full computational power of the 
model. Here we study the behaviour of two non-interacting particles performing a 
quantum walk on the line when the possibility of lattice imperfections, in the form 
of missing links, is considered. We investigate two regimes, statical and dynamical 
percolation, that correspond to different time scales for the imperfections evolution 
with respect to the quantum-walk one. By studying the qualitative behaviour of three 
two-particle quantities for different probabilities of having missing bonds, we argue that 
the chosen symmetry under particle-exchange of the input state strongly affects the 
output of the walk, even in noisy and highly non-ideal regimes. We provide evidence 
against the possibility of gathering information about the walkers indistinguishability 
from the observation of bunching phenomena in the output distribution, in all those 
situations that require a comparison between averaged quantities. Although the spread 
of the walk is not substantially changed by the addition of a second particle, we show 
that the presence of multiple walkers can be beneficial for a procedure to estimate the 
probability of having a broken link. 


I N the field of quantum walks many results have been recently obtained, from both a theoretical [1-5] and an 
experimental [6-14] points of view. Being the quantum analogue of the classical random walk, the quantum walk 
[15] is a simple yet interesting model that can be exploited in several ways. It has been proved that quantum walks 
can be used to implement quantum search algorithms [16], and they can be considered as a universal computational 
primitive [17, 18]. They have applications in the simulation of biological processes [19], and an equivalence with 
quantum lattice gas models [20] has been proposed. It is worth to keep in mind that the scenario with only a 
single walker is however classically simulable, for example using coherent light instead of a single-photon state in an 
optical implementation [21, 22]. In order to have a true quantum behaviour, that could possibly lead to computational 
protocols achieving an exponential speedup over their classical counterparts, the complexity introduced by the presence 
of more than one walker has to be considered (see for instance the short discussion in [23]). From this perspective, 
the study of the emerging features of multi-particle quantum walks, together with an analysis of their resilience to 
non-perfect conditions (e.g. noisy or disordered environments), is of paramount importance. For example, it is well 
known that a quantum particle in a static disordered media cannot evolve arbitrarily far from its starting point, a 
phenomenon known as “Anderson localisation” [24] . Such behaviour has been observed also for particles performing a 
quantum walk, and their localisation length does not considerably change even in the presence of interactions [25-27]. 

In this manuscript we focus on discrete-time quantum walks of non-interacting particles. This model is based on the 
repeated application of a single-step evolution operator, applied independently on both walkers. Such multi-particle 
scenario has been studied when a fixed unitary evolution is considered [21, 28-33]. Since this could not always be the 
case, it is interesting to study what happens with less stringent requirements. In the following we will consider the 
possibility of having missing links in the underlying structure (percolation graphs), modelling for example physical 
situations where one does not have perfect control over the lattice in which the walk takes place. Up to now, this 
scenario has been investigated by several authors only in the single-particle case, approaching the problem from 
different perspectives: e.g. studying decoherence [34] or disordered systems [35], modelling transport phenomena [36], 
or describing asymptotic behaviours [37]. 

The purpose of this paper is to understand how the possible absence of links (due, for instance, to a non-perfect ex¬ 
perimental setup) could influence the evolution of multiple walkers, a fundamental step toward a complete exploitation 
of the quantum-walk model. In order to do so, we study the walk of two non-interacting particles, using as underlying 
lattice a simple line with bonds randomly missing with time. More precisely, on the one hand we want to see to what 
extent signatures of the presence of multiple walkers survive when percolation is considered, while on the other hand 
we are searching for possible advantages arising when dealing with two particles instead of one. To achieve the first 
goal, we qualitatively describe how different quantities depend upon the probability of having missing bonds and on 
the frequency with which the lattice changes. In particular, we will consider the final average distance between the 
walkers and the probability of finding both particles at the same location or exactly in their initial position, showing 
how their values strongly depend on the considered initial state (potentially entangled) even in noisy regimes with 
many broken links. Concerning the advantages of having a second walker we show that, although the spread of the 
walk is not considerably altered by the addition of a particle, there is a procedure to estimate the probability of 
having a missing link whose performance is often improved by its presence. Furthermore, from our study it emerges 
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that a signature of enhanced bunching between the walkers in the output probability distribution does not necessarily 
imply the presence of entangled states symmetric under particle-exchange. Since the symmetry condition is required 
by quantum indistinguishability, this observation is a warning against the possibility of gathering information about 
indistinguishability from the mere presence of bunching effects, at least in all those situations that deal with a random 
component in the evolution by considering averaged quantities. This is an important feature that has to be kept in 
mind also when studying more general many-walker scenarios in future research. 

Let us point out that similar studies in more complex geometries, as higher-dimensional lattices [38] or graphs 
[39], could be of interest. However, being more computationally demanding, this is left as a possibility for future 
investigations. In the remainder of this introduction we will review some basic facts about quantum walks and 
imperfect lattices. 

Single-particle discrete-time quantum walk without percolation. Let us briefly review the concept of discrete¬ 
time quantum walk on the line. This can be thought as the quantum counterpart of a classical random walk, where 
at each step the walker can move left or right with a certain fixed probability. In this scenario, each step can be 
considered as composed of two parts: a (possibly biased) coin toss and a shift in the direction associated with the 
random output. Initially introduced in [15], for a more detailed discussion on the quantum-walk model we refer to 
[ 1 , 2 ]. 

Quantum mechanically the position of the walker is described by orthonormal quantum states that span the 

Hilbert space Hp, while the coin can be taken into account by adding a two-dimensional Hilbert space He, so that 
the global system can be described in Hp <S> He- The coin toss can be realised by means of a fixed unitary operation 
applied to the C component of the system. Throughout this manuscript, we will use for this purpose the Hadamard 
matrix 


H ° y/2 ( 1 “I ) ’ (1) 

expressed in the standard basis of the coin {|t), |i)}, whose states correspond to the classical idea of “head” and 
“tail”. This matrix represents the simplest choice for an unbiased walk, with no direction privileged by the unitary 
evolution. Although it has been shown that the Hadamard coin leads to all possible quantum walks if different starting 
states are considered [40], such equivalence is not guaranteed when percolation is introduced. Our choice has hence to 
be considered as motivated by the sake of simplicity. The second part of a single step corresponds to the shift of the 
walker in the lattice, toward a direction depending on the coin output. Formally, this can be written as the following 
operator acting non-trivially on the whole space Hp <£) He- 

'5 = E(l i + l)< i !®!t)(tl + |i-l)(i|®Wai). (2) 

where we wrote the position component in front of the coin component. We will keep the same order throughout the 
manuscript. 

A single step of the evolution can therefore be obtained by applying the operator 

U = S {tp^Hc) (3) 

to the vector |^ n ) G Hp ® He, representing the initial state of the system. In order to obtain a different result with 
respect to a classical random walk, a measurement is performed in the position degree of freedom only after a certain 
number of steps N. The probability distribution characterising the result will be 


Hi 


,)] (i) = Tr [(|*> (i\ ® 1 C ) UW |VW Win | C/ (JV)t 


(4) 


where the subscript “1” reminds us that it represents the probability associated with a single walker. Without 
percolation, the TV-step evolution operator is simply obtained by exponentiating N times the single-step operator 
U defined in Eq. (3). After the measurement, that instance of quantum walk will be considered as concluded with 
output given by (4). It is well known that the average distance travelled by a classical random walker scales with the 
number N of steps as V^, while for a quantum walker it scales linearly with it [2] . A typical plot for the probability 
distribution of a quantum walk can be observed in Fig. la. 

Notice that all coefficients appearing in the evolution operators (1) and (2) are real. This implies that the real 
and imaginary parts of the initial state (expanded in the basis (|i) |t), \i) \i)}) evolve independently. Therefore, a 
conjugation of the state amplitudes does not change the final probability distribution, a fact that will be used later 
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FIG. 1: Average probability distribution for a quantum walk on the line after N = 300 steps, without missing bonds (a) and 
with dynamical (6) or statical (c) percolation (different ranges for the position values were chosen in order to increase the 
readibility of the plots). In order to have a symmetric distribution we started from the origin with a coin state given by Eq. 
(11). In the last two plots, we considered lattices characterised by a percolation parameter p = 0.75 and we averaged over 10 5 
different outputs. 


Percolation in ID quantum walks. The situation previously described changes if we consider a lattice whose 
bonds can be randomly missing with time. Typically, in a percolation graph each bond is actually present only with 
a fixed probability 0 < p < 1, called “percolation parameter”, and the structure is periodically changed respecting 
such probabilistic constraint. For more information on percolation lattices and their applications we refer to [41, 42]. 
Depending on the frequency of the changes, there can be two extreme regimes of percolation: statical (different lattices 
for different walks) or dynamical (different lattice at each step). The first scenario can be associated with a time scale 
for the lattice-imperfections evolution that is long with respect to the typical time required to finish a single walk, but 
small with respect to the average time needed between different runs of the apparatus. On the other hand, the regime 
of dynamical percolation appears when the imperfections in the lattice evolve so quickly that the missing bonds can 
change position even between consecutive steps of the same walk. In the following both these extreme cases will be 
considered. 

In the presence of percolation, the evolution outlined in Eq. (3) must be changed to avoid crossing of missing bonds. 
To achieve this the step operator applied on the i-th site is chosen as one of the following: 

• S = \i + 1} (i| 0 It) (Tl + I* - 1) (*| ® It) <4-1, 

• S+ = | i) (i\ 0 |t) (t| + |* - 1) (*| ® It) (tl, 

• S- = \i + 1) (i| 0 |t) (tl + N) (*| ® It) (tl, 

• S± = \i) (i\ 0 |t) (tl + K) (*| ® It) (t|, 

respectively when both neighbouring bonds are present, when the following or the previous one is missing, or when 
they are both missing. 

Every time the lattice changes, the missing-bonds positions are chosen randomly, so that eventually the output 
probability distribution has to be averaged over many different lattice sequences. As a result, the ballistic spreading 
typical of a quantum walk is lost, and the final probability distribution approaches that of a classical random walk 
[34, 35], similarly to what happens with other decoherence models [3, 43, 44]. An example can be seen in Fig. 1, where 
the outcomes of a 300-step walk are shown in different percolation regimes. While in the standard quantum walk 
only odd (or even) positions can be occupied, there is no such constraint when the aforementioned strategy of time 
evolution with missing bonds is adopted. Notice also how in the dynamical case the spread of the Gaussian-like shaped 
distribution of the walker is much larger than the one of the statical output probability, which is narrowly peaked. This 
is a typical feature emerging in the two percolation regimes after a certain number of steps given approximately by 
L(p)/ 2, where L(p ) = p/( 1 —p) is the average length of a connected segment of the line, representing the length scale 
over which percolation becomes relevant. After this point the walker evolves with a diffusive behaviour in a dynamical 
percolation regime [34], while in the statical one it remains always localised within a segment: the subsequent average 
leads to the narrowly peaked distribution of Fig. lc, independently of the number of steps. 

Non-interacting particles. In this last introductory section we present the formalism adopted throughout the 
paper to describe two non-interacting walkers. Using the same model of [29], the global Hilbert space of the system 
is l~Li 2 = Hi (8) where each 7~Lk can be expanded as Hk = 'Hp k <S> 1~Lc k - The single-step operator characterising 
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the evolution is chosen to be the single-particle operator U defined in Eq. (3) (or modified close to the missing bonds 
according to the previous rule) applied independently to both systems: 


bi_ L2 = Ui U 2 , 


( 5 ) 


where the label represents the Hilbert space of application. This choice guarantees that both particles evolve inde¬ 
pendently, without any kind of interaction. If the evolution takes place on a perfect lattice, is obtained from 

(5) via exponentiation, i.e., u[^ = U This changes in the presence of dynamical percolation, because every step is 
characterised by a different step operator (associated with a different lattice). 

In particular, in this manuscript we want to focus on indistinguishable particles, satisfying boson-like or fermion-like 
statistics or just being classically impossible to discriminate when detected at the end of the walk. A comparison of 
the corresponding results will allow us to distinguish the features due to the quantum nature of the particles from 
those coming from the averaging process. Two single-particle states \^i) and ^ 2 ) can be combined to describe two 
bosons or two fermions by taking as joint initial state their opportune symmetrisation: 


Y 2 (l ±\(ip 1 \ V>2> I 2 ) 


(6) 


where the indeterminate form |Sym_(^, '0)) 12 is set to 0 by definition. The case of classically indistinguishable 
particles, instead, can be modelled by considering an initial separable state \'ipi} 1 <S> ^ 2)2 an d by changing the mea¬ 
surement stage so as to make the final probability symmetric under particle-exchange. This is done by considering a 
measurement involving the set of projectors 


H-ij ~ 


— Ilji — ^ (\ij) (ij\ Pli p 2 + \ji) 0^Ipi,p 2 ) 


that leads to a probability distribution 


V 2 [IV’in)] (i,j) = Tr \llij l4P | 'Ipin) 


( 7 ) 


(8) 


Notice that when dealing with bosons or fermions the addition of the symmetrised term in the projector does not 
affect at all the output probability. With definition (8), the probability of finding one particle in position i and the 
other one in j ^ i is given by 2 V 2 [IV’m)] (i, j), while the probability of finding them at the same location i is just 
V 2 [IV’m)] (m)- Such description is chosen in order to obtain an output distribution properly normalised to 1 on the 
whole plane (i, j). Given two single-particle states |the output probability distributions obtained from (8) in 
the three cases (bosons, fermions or classically indistinguishable particles) can be labelled respectively by and 
V^ l \ The dependence upon the states l^/c), k = 1,2, will be dropped most of the time to simplify the notation, 
when there is no possibility of confusion. While the first two distributions cannot generally be expressed in terms of 
single-particle outputs (4), for two classically indistinguishable particles the following equality holds: 


j) = ZPi [IV>1>] (i)P 1 


1 


b 2)\ (j) + \pi [IV2)] (i)Pi m] ( j ), 


( 9 ) 


as can be easily proved by using the factorised structure of both evolution (5) and projector (7). 

It is worth noticing that not every entangled state in H 12 can be expressed in a symmetrised form as in (6). An 
example of bipartite state left out by that particular structure is 


^(NotSym)^ _ 


I ij) + \ji) _ in) + in) 


V2 


V2 


( 10 ) 


when i ^ j. In the following, however, we will restrict our investigation to states of a symmetrised form as in (6), or 
to separable states |^i) ® IV^), because we need a definition of k = 1,2, to compare with p!f l \ Moreover, 
in order to maximise the multi-particle effects in our interaction-free model, we will always consider states with a 
localised common initial position: the origin. In addition to this, to avoid asymmetrical spreading, we will focus 
mainly on combinations of the two single-particle coin states 


\<p±) = 


It) ± i It) 

V2 ' 


( 11 ) 










5 



-10 


-15 


FIG. 2: Output probabilities for a 15-step quantum walk of two bosons, fermions or classically indistinguishable particles 
[respectively in (a), ( b ) and (c)], starting from the origin with the coin states given in Eq. (12). 


that lead to the symmetric evolution shown in Fig. 1. With these conditions, the corresponding coin states for bosons, 
fermions and classically indistinguishable particles read 


\<f>+) = 


in) + in) 

V2 


in) 


in) - in) 

V2 


in) = w+) ® W-)- 


( 12 ) 


Notice that, with the initial state IVts'), relation (9) further simplifies: being \tp+) and \(p~) in Eq. (11) the complex 
conjugate of each other, their single-particle output distributions are the same: 


V 1 [\v+)}=1-\[W-)]=V l , (13) 

so that the joint probability becomes 

v[ cl \i,j)=V 1 {i)V 1 {j). (14) 

Examples of the probability distributions obtained by using the initial states defined in Eq. (12) can be observed in 
Fig. 2, where typical bunching and anti-bunching behaviours can be noticed (see [33] for an experimental realisation). 


Results 

We can now finally address the problem of considering the quantum walk of two indistinguishable and non-interacting 
particles in a one-dimensional percolation lattice, where every run of the walk (in the statical case) or even every 
step (in a dynamical regime) takes place on a different underlying structure, randomly generated as described above. 
Considering a number A E N of quantum walk realisations, the average output probability distribution reads 

_ ^ 

(*’•?)’ ( 15 ) 

a= 1 

where index a runs from 1 to A. The evaluation of such quantity for different input states will be the basic tool 
adopted to characterise and study the behaviour of the walk. Usually the plots in this manuscript are obtained by 
averaging over a number of different outputs varying from A = 2000 to A = 5000, depending on the precision required 
and on the complexity of the simulation. The obtained results do not qualitatively change in a significant way in that 
range. Moreover, they do not strongly depend upon the number of steps considered, that will typically be N = 15 
(we will explicitly point out the features that depend upon the parity of N). Therefore, we can use the obtained 
information to draw conclusions about the general qualitative behaviour of the walk. Our discussion will include three 
sections, addressing three different issues. 

At first, we want to qualitatively describe the average output distribution obtained in different percolation regimes 
considering the symmetric input states of Eq. (12). The goal of this analysis is to gain a better intuitive understanding 
of the situation, pointing out the features that are due to the quantum interference between the particles (that arises 
for entangled initial coin states, associated with bosonic or fermionic statistics) and those effects whose presence 
is rooted in the average appearing in Eq. (15). In particular, we show how a signature of enhanced bunching 
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FIG. 3: V\ (i,j) for a quantum walk of two walkers starting from the same position with coin state | ips) [see Eq. (12)], 
measured after 15 steps. Panels (a) and (5) correspond respectively to a dynamical and statical percolation regime, with 
p = 0.9 and a number of averages A = 1000. 


between the walkers in the output distribution does not necessarily imply the presence of entangled states symmetric 
under particle-exchange. This is therefore a warning against the possibility of gathering information about quantum 
indistinguishability from the mere presence of bunching effects. 

In the second part of our discussion, we study how some quantities of interest depend upon the parameters of the 
problem (percolation regimes, probability of missing bonds and input states). In particular we focus on the average 
distance between the positions in which the particles are detected (D) [29] and on the probability of finding them 
either in the same location (M), also called “meeting probability” [30, 45], or exactly in the origin (C). On the one 
hand, this aims at capturing some particular aspects of the output by means of a few parameters that can be more 
easily interpreted and manipulated with respect to the whole probability distribution. Being the whole distribution 
impossible to plot for higher dimensions of the lattice or for a larger number of particles, in future research this kind of 
approach will become more and more common. From this perspective, our analysis aims at highlighting the qualitative 
trends that appear in the aforementioned quantities in different scenarios of interest, whose identification will ease the 
interpretation of future, more complex, results. On the other hand, we are interested to see if the behaviour of these 
quantities strongly depends upon the considered input state [chosen among (12)], implying the relevance that inter¬ 
particle quantum interference maintains even in very noisy regimes. This fact, observed in particular for the meeting 
probability M, can be interpreted as a signature of resilience of multi-particle effects to structural imperfections. 

Finally, in the third and last section, we search for possible advantages arising when considering a multi-particle 
quantum walk with respect to its single-particle counterpart. At first we consider the single-particle spread, showing 
that it cannot be improved by the addition of another particle. However, inspired by the qualitative behaviour of 
C that was previously found, we show how such quantity can be used in a procedure to estimate the percolation 
parameter (probability for each link to be present) whose performance can be improved by the presence of a second 
walker. 

_ ( c /) 

A. Analysis of bunching events. The typical behaviour of the average output distribution P 2 for classically 
indistinguishable particles, defined by evaluating Eq. (15) with the separable input state | ips) of Eq. (12), is plotted 
in Fig. 3 for dynamical and statical percolation. As can be appreciated, the two results are qualitatively dissimilar. 
In particular, we can notice how the distribution is more spread out over the allowed region in the dynamical case, 
while in the statical regime it is concentrated in a few particular points, mainly along the diagonal. This result can 
be interpreted by keeping in mind that in the first case the walkers spread in a diffusive way, while in every run of 
a statical regime they are localised on a small segment around the origin. The final distribution is then obtained by 
averaging the output of each run, whose detection probability is similar to Fig. 2c. This leads to the diagonal peaks 
visible in Fig. 3b, where the corresponding off-diagonal peaks of Fig. 2c disappeared because it is very unlikely to 
have a connected segment that extends far away from the origin in both directions. Similar considerations could be 
made for entangled coin states, representing boson-like and fermionic-like statistics, if instead of Fig. 2c we consider 
as underlying coherent evolutions for example those given in Fig.s 2a and 2b. From this analysis it emerges that 
the so called “bunching events”, associated with the detection of two particles at the same location, depend on the 
interplay between statistical average and effects of quantum interference between multiple walkers. In the following 
we want to formalise these considerations. 

In order to properly define a peak (or a valley) in an output distribution, we need a reference probability to use 
as comparison. The peaks (or valleys) due to inter-particle quantum interference can be evidenced by comparing the 
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average probabilities V^\ associated with the symmetrisation of two single-particle input states |-0i) and 02 ) [see 

Eq. (6)], with the average probability V 2 , corresponding to classically indistinguishable particles starting from a 
separable input state obtained by taking the tensor product of the same \^\) and IV^)- Indeed, for a single run of the 
quantum walk, whose evolution is labelled by a, we can expand the probability of finding both particles in the same 
location as 


v ( 2 ±Ka \j,j) 


1 ± | (V>i| 1P2) I 2 


± 


(V’i 


li)p 


01 u{a) IO 2 ) 


1 ± I (V’iI ^2) I 2 


(16) 


From this, it follows that at least when (Vq| ^ 2 ) = 0 [which is the case in the main situation of interest here, where 
|^ 1 , 2 ) = |0> |<p±)] the elements on the diagonal are larger (for bosons) or smaller (for fermions) than the diagonal 
elements in the corresponding probability for classically indistinguishable particles, as can be expected. Since this 
property holds for every index a, it is still true for the average probability distribution. This effect in a regime of 
dynamical percolation can be observed in Fig. 4. 

_ (cl) 

On the other hand, the effect of the average can be understood by comparing V 2 with the symmetrised product of 
the two single-particle average distributions V\ [|V^)]- Surprisingly, in some situations the averaging procedure might 
induce an enhancement of the probability of bunching events, leaving a signature in the output distribution similar 
to the interference effect previously analysed. To see this, we can expand the probability v[ a ^ [\^k)] around its mean 
value as 


v i c ) m\ = v,m]+5^u k )], a?) 

so that, for every lattice site j, the differences 5^ [I'lpk}] average to zero: 

jJ2s^[\^ k )](j)=0, ( 18 ) 

< 2=1 

giving, for the average detection probability of two classically indistinguishable particles, 

pf\hj) = \{p tdV’i))] w p [(i^2»] a) + p KiV’i))] u) p im)} (*)) 

+ ^ E ( s(a) [( 1 ^ 1 ))] (o -5 (a) [d^2»] u) +<5 (a) KiV’i))] u) s (a) im)] (*)). ( 19 ) 

a=l 

Usually when \ipi) and ^ 2 } lead to different probability distributions, or when off-diagonal terms i ^ j are considered, 
the second line in (19) has not a fixed sign, so for A 1 it does not contribute significantly. This is no longer true 
when we consider diagonal terms i = j for states with the same final distribution (e.g. |0) |</?+) and |0) |<^_)). In this 
case the second line in (19) is always positive, being a sum of squared terms. These contributions hence add up to 
create the peaks over the first term, which in the example provided is just equal to the product (see Eq. 

(14) and Fig. 5 for the resulting plots in a regime of dynamical percolation). 

We expect that a similar behaviour can arise also with more complicated lattice structures affected by percolation. 
Therefore, whenever one deals statistically with such imperfections, it is useful to keep in mind that the mere av¬ 
eraging process could lead to an enhancement of the diagonal entries of the average probability distribution. Being 
aware of this fact can avoid the misinterpretation of such phenomenon as a bunching signature due to the quantum 
indistinguishability of the multiple walkers considered. 

B. Qualitative analysis of three significant quantities in the presence of percolation. In what follows, 
we want to see how a few quantities depend upon the adopted percolation parameter p and the particle-exchange 
symmetry encoded in the input states, for both dynamical and statical percolation regimes. In particular, we will 
focus on two walkers starting from the same position with initial coin states | </>+), IV’-) or | defined in Eq. (12). 
For every considered quantity Q and for a set of equally spaced percolation parameters p, we will average the value 
obtained for a single lattice realisation Q^ over A = 5000 different configurations: 

1 A 


( 20 ) 
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(d) Difference (a) — (c) (e) Difference (6) — (c) 


FIG. 4: Average output probability distributions for a two-particle 15-step quantum walk affected by dynamical percolation 
(p = 0.75, A = 2000). We considered boson-like (a), fermion-like (5), and classically indistinguishable (c) combinations of two 
single-particle states \p±) that spread symmetrically from the origin. Panels (d) and (e) evidence the contribution to bunching 
events due to inter-particle quantum interference. 




(b) Pi(0 Pit i) 



p 


(c) Difference (a) — ( b ) 


_( c n _ _ 

FIG. 5: ' (h j), (j), and their difference are plotted respectively in (a), (5), and (c), for a quantum walk of N = 15 

steps starting from the origin with an initial coin state {'ips ), defined in Eq. (12). A dynamical percolation regime with p m 0.75 
averaged over A = 2000 different evolutions is considered. This figure shows how the statistical averaging procedure might 
enhance bunching events in the average probability distribution of a multi-particle quantum walk. 


The uncertainty associated with such average will be then given by the standard error 


TQ (P) = A^l ^ ~ G(P)f 


( 21 ) 


Bl. Final average distance between the particles. A global property that can be analysed is the average 
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FIG. 6: Average output distance between the two particles as a function of the percolation parameter p, for a 15-step quantum 
walk starting from the origin with three different coin states [see Eq. (12)]. (a) and ( b ) correspond respectively to dynamical 
and statical percolation regimes. The average is calculated over A = 5000 different lattice realisations and the errors are 
evaluated as the standard errors on the mean values obtained. 


distance between the positions where the two walkers are detected. Formally such quantity can be expressed as 

D = Y^\3-i\V 2 (i,j). ( 22 ) 

*3 

Its numerical evaluation for different percolation parameters p is plotted in Fig. 6, for a 15-step quantum walk. As 
can be expected, for any value of the percolation parameter, this quantifier for fermions is larger than the one for 
bosons, with an intermediate behaviour for the separable coin I'lps}- This result agrees with the discussion in [29], 
that considered this same quantity without percolation. 

If we consider the dependence upon p in the dynamical percolation regime, we can see that the average distance 
between bosons increases approximately in a linear way, with a small decrease just before p = 1. This can be 
interpreted as a trade-off between two trends: on the one hand with less missing bonds the two walkers can move 
on a larger segment of the line, increasing the likelihood of finding them further away; on the other hand, if the two 
particles are less disturbed, their distribution approaches the one obtained without percolation (shown in Fig. 2a), 
characterised by a high degree of bunching. For fermions, instead, the same considerations add up increasing the 
average distance of separation, since they tend to show an anti-bunched distribution in the limit p 1 (Fig. 2b). 

If a comparison is made between the plots obtained with dynamical and statical percolation for a 15-step walk, 
one can notice that in the second regime the average distance between the particles is much smaller. This can be 
understood by thinking that with statical percolation in each run the particles are constrained to stay close to each 
other in a fixed region of space, whereas in the dynamical case they have always the chance of getting further away. 
The observed behaviour is hence typical once the number of steps is larger than the average length of the connected 
segment L available to the walkers in a regime of statical percolation [34]. 

B2. Meeting probability. Another global quantity, easier to obtain with respect to the previous one as it involves 
a smaller number of measurements, is the overall probability for the two walkers to be on the same position at the 
end of the walk. Such quantity, also called meeting probability [30, 45], corresponds to the sum of the diagonal terms 
in the output probability distribution V 2 (i,j) in Eq. (8): 

M = (23) 

3 

The results obtained for a 15-step quantum walk are plotted in Fig. 7 for both dynamical and statical percolation. 
We can notice how in the first scenario the meeting probability decreases for p < 0.92 and slightly increases after 
this point, while in the statical case it is monotonically decreasing for all values of p. The minimum position in the 
dynamical case seems to depend on the number of steps, moving slowly toward p — 1 when this is increased. However, 
even if this could look like a finite-size effect, we pushed simulations up to 80 steps, and have strong numerical 
evidence that the non-monotonic behaviour would still be present even for a much larger number of steps. We can 
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P 


(b) 


FIG. 7: Average meeting probability as a function of the percolation parameter p, for a 15-step quantum walk starting from the 
origin with three different coin states [see Eq. (12)]. (a) and ( b ) correspond respectively to dynamical and statical percolation 
regimes. The average is calculated over A = 5000 different lattice realisations and the errors are evaluated as the standard 
errors on the mean values obtained. 


obtain an intuitive explanation for this fact by comparing the output of the walk in a perfect lattice (Fig. 2) with the 
“spreading” effect that dynamical percolation has on the output probability distribution (see Fig. 3a). Starting from 
the ideal case of p = 1, the addition of small amounts of imperfections makes the diagonal peaks of Fig. 2 spread, 
leading to a diminished meeting probability. This effect is much less important in a statical percolation regime, as 
can be appreciated in Fig. 3b. This is because in every lattice realisation the walk is coherent on a small segment of 
the line around the origin, and the percolation only averages such small-scale coherent evolutions. The overall result 
is that in this case the meeting probability depends monotonically upon p (see Fig. 7b). 

From Fig. 7 we can also see how the quantum statistics of the particles heavily affects the meeting probability, 
even for small values of p. This shows that a signature of bosonic or fermionic multi-particle effects survives also in a 
noisy regime where the quantum walk tends to its classical counterpart. 

B3. Probability of finding both walkers in the origin. Finally, we consider a “local” property, i.e. that can be 
obtained through a measurement in a single position: the probability of finding both walkers in their initial starting 
point. In the following we will label it as 


C = v 2 { 0 , 0 ). (24) 

In addition to being more easily obtained with respect to the previous quantities, this one differs also because it 
depends on the parity of the number of steps N considered. Indeed, when p = 1, C is rigorously zero if N is odd, 
but this could not be the case if N is even. Being mostly interested in this quantity for its monotonic and regular 
behaviour, that can be exploited to estimate the percolation parameter p (see next section), we will focus only on odd 
numbers of steps in order to have a quantity varying in the maximum range [0,1]. 

Intuitively, C will be close to 1 for small values of p, since both walkers are constrained by the missing bonds to 
stay close to the starting point, and it decreases with increasing p while the walk spreads further away, reaching zero 
for p = 1 (for N odd). This behaviour is numerically confirmed by Fig. 8, that shows how C depends upon p for a 
15-step quantum walk affected by dynamical or statical percolation. We can see how in the first regime C(p) goes 
quickly to zero as soon as the imperfections are reduced. This can be explained with the diffusive behaviour of the 
walkers, that are allowed to quickly populate regions with i ^ j in the V 2 (h j) plot (e.g. see Fig. 3a). The situation 
is different in the statical case, where the walkers move away from the origin in a more controlled way (e.g. see Fig. 
3b), based once again on the typical length of connected lattice around the origin. The resulting plot for C(p) is 
monotonic with almost constant slope in the whole range [0,1], as can be appreciated in Fig. 8b. 

C. Search for multi-particle advantages. In this final section we want to discuss whether the addition of a second 
walker to the walk can lead to an immediate advantage in some task or protocol. As a first check we consider the 
possibility of having an improvement in the spread of the walk and show that this is not the case. For example, 
for every quantum statistics of the two particles involved there is always a single-particle state leading to a larger 
spread. We then consider a more specific task, namely the estimation of the percolation parameter affecting the lattice, 
proposing a two-walker procedure that in certain regimes shows an advantage over its single-particle counterpart. 
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C(P) 



(a) 


C(P) 



FIG. 8: Average probability of finding both walkers in the origin, as a function of the percolation parameter p, for a 15-step 
quantum walk starting from the origin with three different coin states [see Eq. (12)]. (a) and (5) correspond respectively to 
dynamical and statical percolation regimes. The average is calculated over A = 5000 different lattice realisations and the errors 
are evaluated as the standard errors on the mean values obtained. 


Cl. Spread of the walk. The intuitive idea of spreading for a single walker, starting from a localised position zq, 
can be well described by the quantity 

Vi = ^(i-io) 2 J Pi(i), (25) 

i 

where the initial coin state is not explicitly written. Notice that, if the walker evolves symmetrically, such quantity 
corresponds to the variance of the output distribution. How can we generalise this expression to the two-walker case? 
Also here, we are mostly interested in letting both particles start from the same lattice position. Indeed, since in 
our model there are no interactions, this situation is the one that could maximise the multi-particle effects. For this 
reason, we will consider here only the case where both walkers start from the origin (with arbitrary coin states). With 
this hypothesis, the straightforward generalisation of (25) for a generic fixed unitary evolution labelled by a (i.e. for 
a given sequence of lattices) is 

F Q) = ^E( i2 +G^2 a Ai)- (26) 

ij 


Once two single-particle coin states have been fixed, we can label with the spread corresponding to 

the different statistics of the particles, obtained by using in (26) respectively or p( cl ^ a \ In the same way, 

with an overbar we will label the value of the spread averaged over many different lattice sequences in a percolation 
regime. The overall factor 1/2 in our definition (26) is introduced only to ease the comparison with the single-walker 
spread V\ of Eq. (25), because from Eq. (27) it follows 


(cl)(a) _ F 0) (Cl)+VT( c 2 ) 




(a)( 


< max Vi (a) (cfc), 

ke{ 1 , 2 } 


(27) 


where two generic single-particle coin states |ci), \ci) are considered. This is a simple consequence of the fact that 
two particles initialised in a separable state spread independently under the evolution in Eq. (5). The fact that we 
cannot distinguish them in the measurement leads only to the average that appears in Eq. (27). 

One can now wonder whether the presence of a proper quantum statistics, introduced by the symmetrisation of the 
vector describing the state, can somehow alter this result. In order to show that the situation does not considerably 
change, we consider an initial state 

^in ) ) = l°°> ® |Sym ± (ci,c 2 )> , (28) 


whose spread after N steps (for a particular lattice sequence) can be written as 


k =1 


( 29 ) 








12 


In this expression 0 < A 


(±) 

k 


< 1 and 



G Hc 1 are the two eigenvalues and eigenvectors of the coin state 


Pi 


(0|Xr 2 






(±) 


|0>j 


(30) 


obtained by tracing out the Hilbert space of the second particle. Moreover, the odd combination |Sym_(ci, c 2 )) of 
two generic coin states \ci) , |c 2 ) G Tic, as well as the even one |Sym + (ci, C 2 )) or separable one |ci) (8) |C 2 ) of two 
orthogonal coin states (ci| C 2 ) = 0, yields the same spread as the two single-particle symmetrical-spreading states 
| ip±) [defined in Eq. (11)]: 




(-)(«) = y(+)(a) 


<ci| c 2 )=0 


_ pd c 0( a ) 


<ci| c 2 )=0 


= v} a \v±). 


(31) 


The proof of these facts can be found in Methods. Being | ip±) orthogonal coin states leading to the same probability 
distribution, Eq.s (27) and (31) imply 


d a) (</>+) = =v^^s)= y a >±)- 


— T/( a ) 


— T/( a b 


— T/( a b 


(32) 


This means that there is no advantage, from the spread perspective, in considering two symmetrical-spreading particles 
instead of just one, no matter what the statistics over particle-exchange enforced in the input state by Eq. (6) is. 
From Eq. (31) it follows that only the bosonic statistics could a priori influence the spread, when two non-orthogonal 
single-particle coin states \c\) and |c 2 ) are symmetrised. However, Eq. (29) implies that there will always exist a 
single-particle state with a spread at least as large as the one found with two bosons, because 


T/( + )( a ) 

v 2 


Zx?>v<‘>(t , )< t ™ S) v 1 < " ) ( 4 +> ) ■ 


(33) 


We emphasise that such states 


C ( , +) ) 


w k I could be different from both coin states \ci) and \c 2 ) symmetrised at the 

beginning, and that the optimal k could depend on the explicit lattice realisation labelled by “a”. 

All these results can be trivially extended to the spread averaged over many lattice configurations. In particular, 
for every initial coin states \ci) , \c 2 ) one has 


jr(cl) Vi(ci) + V\ (c 2 ) 


n = 


< max Hi (cfc), 
ke{ 1 , 2 } 


W = E 4 + E (4 +) ) <_ max. Vi 


k~l 


ke{ 1 , 2 } 


( 4 + >), 


^2 } =V 1 (<p ± ), 


in terms of 


4 +) ) and A^, + ^ defined in Eq. (30). 


(34) 


(35) 


C2. Estimation of the percolation parameter. Here we want to explicitly consider the possibility of using 
some of the quantities described in section B to infer the value of the percolation parameter p that characterises the 
lattice imperfections. In order to do so, let us consider a certain event E, happening at the measurement stage of the 
quantum walk: e.g. both particles at the same location, both in the origin, etc. The probability of such event will 
depend on the explicit lattice realisation, and we will write it as PgP, where £ runs over all the possible sequences of 
lattice configurations. Notice that here we are using £ instead of the previous a, to emphasise the fact that we are 
considering an effective lattice realisation and not of a simulation of it. has to be interpreted as the probability of 

the event considered, conditioned upon the £-th realisation of the lattice: = P(P|£). Therefore, the probability 

of E can be written as 

Pe = E P(P\0 P(0 = E p{ E (36) 

£ £ 

This quantity, although difficult to obtain analytically, can be numerically estimated. Indeed, for the event “both 
particles in the same place” (or “both particles in the origin”) it would correspond to M(p) (20) [or C(p) (24)] in the 
limit of many different sequences of lattice configurations considered (A —>■ oc). In the following we will generically 
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label this numerical estimate as P^ zrn \ whose precision will depend on the number A adopted [see the standard 
error in Eq. (21)]. On the other hand, if one had access to an experimental setup of such quantum walk affected by 
percolation, after n runs of the experiment Pe could be estimated via the quantity 


p ( est ) _ 

TP - l 


(37) 


where tie is the number of times in which E happened. This estimator is a stochastic variable, characterised by 
expectation value and variance given by 


EV 


p(est) 

r E 


= Pe, 


VAR 


p(est) 

r E 


Pe(1~Pe ) 

n 


(38) 


This can be shown by considering P^ st ^ as the average of n independent and identically distributed variables, taking 
the value 1 with probability Pe (when E happens in that single run) and the value 0 otherwise. 

Our goal is to estimate the value of p through a measurement of Pg St \ that estimates Pe with an error 


VAR 


0 (est) 

E 


1/2 


Exploiting the linear propagation of errors and the explicit expression of the variance given in 


Eq. (38), we can write the uncertainty in the estimation of p as 


/V,AR 

p(est) 

r E 


dP E 



dp 


\APe(1 - Pe) 



(39) 


so that an upper bound Sp < e yields a lower bound on the required number of runs n > r// n (p ): 


\ (e) 

n > n min 


(p) = 


Pe( 1 - Pe) 


dP E 

dp 


(40) 


The right-hand side can be estimated by using P ^ trn ^, numerically obtained by considering a total number of A 1 
sequences of lattices. In particular we can approximate 


Pe( 1 - Pe) - P { E m \ 1 - P { E m) ), 
dP E _ dP ( E sim) 
dp dp 


(41) 

(42) 


where the second term is obtained from a polynomial fit of the numerical samples {P^ im \pi)}i. We point out that, 
despite Eq. (38), the expectation value of the right-hand side of (41) coincides with the left-hand side only in the 
limit A oo, making P^ zrn \l — P^ irn ^) a biased but consistent estimator for Pe(1 — Pe)- 

In order to explicitly obtain the bound n^ in (p) for different kinds of quantum particles, we considered the event 
“both walkers in the origin” for a 7-step quantum walk in a regime of statical percolation. Indeed, having an almost 
constant slope over the whole region p G [0,1], C(p) is the best quantity among those studied in subsection B to 
estimate the value of p in that regime. Notice that here the number of steps N considered is reduced (with respect to 
the usual 15) to allow us to take a much larger value of A: for this analysis we randomly generated ^ 10 5 percolation 
lattices. This guarantees that the total number of simulated configurations A is much larger than the total number 
of possible lattices, given by 2 7 ' 2 ~ 1.6 • 10 4 for the situation considered. We stress once again that in order to obtain 
plots qualitatively similar to Fig. 8b, leading to a good estimation strategy, it is important to consider an odd number 
of steps. 

The single-walker counterpart of the event here considered corresponds to detecting the particle in the origin at the 
end of the walk, and it can be used to infer p with the same procedure previously outlined. Therefore, we can compare 
the performances achieved using either one or two walkers initialised with different input exchange symmetries. In 
Fig. 9 we plotted the event probability P^ irn ^ and the bound for the input coin states defined in Eq. (12), 

corresponding to boson-like, fermion-like statistics or to two classically indistinguishable particles. We can see how 
for the 7-step walk considered, for values of the percolation parameter p < 0.82, there is an advantage in using two 
quantum particles instead of one. In particular bosons are optimal in the range 0.38 < p < 0.82, while fermions 
perform better below p ~ 0.38. We studied the dependence of such thresholds of bosons optimality pi < p < P 2 upon 
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FIG. 9: (a) Simulated event probability P/f Z7T ^ of finding the single particle used (SP) or both particles [see Eq. (12)] in the 
origin, as a function of the percolation parameter p. Due to the large number of simulations, the standard error is negligible. ( b ) 
Corresponding lower bound on the number of required experimental runs as a function of the percolation parameter 

p. In both cases a 7-step quantum walk has been considered, averaging over A ~ 10 5 lattice sequences. For the evaluation of 
n mln> we used polynomials of degree 5 to fit the points in (a) (see the discussion below Eq. (42) in the main text). 


N, considering walks with an odd number of steps between N = 3 and TV = 11. We found that p\ (P 2 ) decreases 
(increases) slightly with increasing AT, with pi and P 2 varying respectively in the ranges [0.35,0.44] and [0.78,0.83]. 
This suggests that the range of bosons optimality could increase with the number of steps considered. 


Finally, let us comment on the limit n^ in —>> 0 that Eq. (40) yields for Pe —» 0 or Pe 1. The error VAR 


A est ) 


1 1/2 


(and therefore Sp) is associated with a certain probability of finding the real value in that range, that for a small 
number of measurements n is actually unknown. Only for n 1 the central limit theorem applies, allowing us 
to interpret the obtained error as based on a well known Gaussian probability distribution. Moreover, the linear 


propagation of errors assumes VAR 


A est ) 


1/2 


small enough for the function Pe{p ) to be considered approximately 


linear in that interval. Even if Eq. (40) seems to suggest otherwise, the scaling ^ 1/ y/n in Eq. (38) therefore requires 
to take n^- large enough to satisfy the previous conditions. 


Discussion 

We studied the behaviour of a one-dimensional quantum walk when the possibility of having missing bonds is consid¬ 
ered together with the presence of a second walker. We focused on two percolation regimes, statical and dynamical, 
associated respectively with slowly and fast varying imperfections. Our numerical simulations have been averaged over 
many (typically ^ 5000 for a 15-step quantum walk) different lattice sequences, realised by associating to each bond a 
probability p < 1 of being present. Starting from two symmetrical-spreading single-walker states, we combined them 
in three different ways, in order to observe the behaviour of two bosons, fermions (showing inter-particle interference), 
or simply two classically indistinguishable particles (with a non-symmetrised quantum state). 

We described how the bosonic (fermionic) quantum statistics of the walkers affects the final probability distribution, 
inducing the presence of characteristic peaks (valleys) in its diagonal values. This feature, easily observed through a 
comparison with the output of an unsymmetrised state, is not due to percolation but emerges as a consequence of 
the imposed exchange symmetry between two orthogonal single-particle states. On the contrary, the remaining peaks 
that may still be observed in the distribution of two classically indistinguishable particles are due to a statistical effect 
of the average over different unitary evolutions, associated with diverse realisations of the underlying lattice structure. 
This shows how in the presence of percolation high probabilities of finding both walkers at the same position do not 
always imply a bunching phenomenon due to the bosonic nature of the particles. 

The advantage of studying low-dimensional lattices with no more than two walkers is the possibility of representing 
the output detection probability in a 3D plot. This visual approach, however, cannot be extended to more compli¬ 
cated scenarios, involving multiple walkers and/or more complex structures. The search for quantities that could 
meaningfully describe some aspects of the walk behaviour, being at the same time easily accessible theoretically and 
experimentally, will therefore be of extreme importance in the future. Here we considered three of these quantities 
(final average distance, probability of finding both walkers at the same location or exactly in the origin) describing 
their qualitative behaviour for a symmetrical-spreading two-particle quantum walk with percolation. In particular we 
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saw how they often depend strongly on the quantum statistics of the walkers even for low values of the percolation 
parameter p, fact that suggests the resilience of inter-particle interference to very noisy conditions. 

Searching for multi-walkers advantages, we considered at first the spread of the walk, focusing on the dependence 
upon the second particle more than on the decoherence introduced by the percolation. Our analysis shows how the 
addition of a second walker does not increase the spread, since there always exists a single-particle state with a spread 
at least as large as the one obtained with two walkers. Moreover, for many cases of interest, i.e. for any couple of 
fermions and for bosons or classically indistinguishable particles obtained by combining two orthogonal single-particle 
coin states, the spread is not affected by the quantum statistics of the walkers. In these cases it is equal to the spread 
of the two symmetric single-particle coin states |<p±) [defined in Eq. (11)]. 

We also proposed a procedure to estimate the value of the percolation parameter that characterises the amount 
of lattice imperfections. Due to the monotonicity and the regularity observed in Fig. 8b, it is worth to apply such 
strategy considering the detection of both particles in the origin in a regime of slowly varying imperfections, i.e., of 
statical percolation. For values of p not too large, we showed how the performance of the proposed estimation is 
enhanced by the presence of a second walker, that allows us to need a smaller number of experimental measurements 
in order to achieve a certain estimation precision. 

This investigation considered for the first time the effects of possible missing links (that could happen, for instance, 
in an experimental setting) on the features of walks where the quantum nature of the involved particles is relevant. 
Due to the importance of quantum walks for the purposes of quantum computation and simulation, and in particular 
for the role that many-walkers complexity will have in the future, these results pave the way for a more complete 
understanding and a full exploitation of these models. 


Methods 

Comparison between single-particle and two-particle spread. Here we will provide a proof for Eq.s (29) and 
(31). In the following the term single-particle spread will always correspond to the quantity V\ obtained from Eq. 
(25) by setting io = 0 and we will use the shorthand notation V% and Vij to represent respectively the single-particle 
and the two-particle output probabilities to detect the walkers in positions i and j. We will also drop the dependence 
upon the lattice realisation “a”. 

Let us begin by showing that given two orthogonal coin states \c\) , fa) £ He , characterising the initial conditions 
of two single walkers localised in the origin, their average spread depends only upon the unitary evolution and 

not on the basis {\ci }, fa)} of He- Indeed, being 

fa )c fa\ + fa )c fa\ = H-Cj (43) 

from Eq. (4) the sum Pi[fa)] + Pi[fa)] does not depend on the choice of basis: 


Pi[\c 1 )]+P i [\c 2 )]=Tr 


C/ (JV)t | i)p (i\ ® 1 c u {N) \i 0 )p Col ® 1 c 


(44) 


property that holds also for the average between the two single-particle spreads. Therefore, given two particles starting 
from the origin with two orthonormal coin states fa ), |ci) £ He'- 


HW + v, W - yi(w ) 


(45) 


being (<p+| <p_) = 0 and V\ (<p+) = Vi(<p_). 

We are now able to prove Eq. (29). Being Vij symmetric under the exchange 2 £ j, the two-particle spread V 2 
defined in Eq. (26) only depends upon the marginal sum JT Vij: 


b ±} = E* 2 E^ 


(±) 
ij ' 


(46) 


Using the symmetry under particle exchange of the evolved state 


^)=u[P 


lb) 

T V 




(47) 


from Eq. (8) it follows 
1 
2 


E P t ±) = 5 E ^ 1 ’ 2 (!*•?’> (*•?! + U*> 0'*l) V’/ ±} ) (V’/ ±) = Tri (|*)i <*| • u[ 


(N) 


Tro 




77^} (48) 


= E4 ) 1*i |i> 1 <i|-^r ) |0) 1 (0| 


km 1 


: (±)\/-(±) 


Vi™ =Z4 ±> r. 


k=l 


;(±) 


(49) 
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where the last two equalities follow from the definition of and ^ as eigenvalues and eigenvectors of (30), 
using 


Tr 2 


4n ] ) (4n\] = |0) Pl <0|®E4 ±} |4 ±} ) Ci (c| 


(±) 

7c 


(50) 


The required equality in Eq. (29) then follows trivially from the spread definition (46). 

Eventually, it can be easily proved that for every \c \), |c 2 ) E Pic one has s ~, k 1, 2, and that the same 
equality holds also for bosons if (ci| c 2 ) = 0. Therefore, Eq. (31) is a consequence of Eq.s (45) and (27). 
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